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Abstract 

Mathematical equations are derived by UBing tlio 
Maximum Prlnolplo to obtain tlie maximum payload ca- 
pability of a reusable Tug for planotary missions. 11 10 
mathematical formulation includes correction for nodal 
precession of Llio Space Shuttle orbit. The Tug performs 
this nodal correollon In reluming to this processed orbit. 
The sample case analyzed represents an Inner planet 
mission as defined by tho declination (fixed) and right 
ascension of the outgoing asymptote and the mission on- 
orgy. Payload capability Is dorlved for a typical cryo- 
genic Tug and the sample cnBo with mid without perlgeo 
propulsion. Optimal trajectory profiles and some im- 
portant orbital oloments aro also discussed. 

Introduction 

At tho present time NASA is developing a rousablo 
Space Shuttle (SS) that will carry pa'-’aads to low Earth 
orbit. For missions beyond low Earth orbit the SS will 
ourry In Its cargo bay (in addition to tho payload) a pro- 
pulsive stago that will be deployed In low Earth orbit and 
will deliver tho payload to Its required injection condi- 
tions. 

This propulsive stage may bo oxponded after pay- 
load delivery, or it may be retumad to a woittng SS for 
a flight back to Earth, Expandable stage performance 
can bo analyzed with techniques dovoloped for expend- 
able launch vehlclos and will not bo discussed horo'n. 

In this paper, the maximum performance capability 
of a rousablo stage (Tug) is derived for planetary mis- 
sions, Tliis Is a complex mathematical optimization 
problem since tho Tug is beyond Earth-escape energy 
at payload injection and must perform a rolrobum to re- 
turn to an Earth orbit. Since this retrobum occurs far 
from Earth, large volocity losses are encountered, as 
will bo discussed. 

The mathematical optimization problem is formu- 
lated by using the Maximum Principle. State and asso- 
ciated adjoint equations aro numerically integrated to 
detorminc tho instantaneous position and velocity of tho 
Tug. To obtain tho maximum payload capability, varia- 
tional final conditions are derived by using tho Maximum 
Principle, and the two-point boundary value problem is 
solved by using a Newton- Raphson iteration technique. 
The required partial derivatives of final conditions with 



respect to initial conditions aro evaluated by porturblng 
tho initial conditions one at a time, Integrating tho state 
and adjoint equations, and observing tbc changes in final 
conditions, 

Reusable Tug trajectories, including nodal correc- 
tion without porigoo propulsion, have been investigated 
In rcfercnco 1, The analysis presented heroin extends 
the rosults of this reference by optimizing tho Tug total 
trip timo and introducing perigee propulsion. Total trip 
timo is defined as tho elapsed timo from tho start of tho 
first outbound bum to return to tho SS orbit for rendez- 
vous. For porlgcc propulsion the continuous outbound 
bum required to roach payload Injection conditions Is 
split into two bums separated by an optimum coast timo. 
Tho present analysis also Includes nodtl regression of 
the SS orbit, Nodal regression Is caused by tho Earth's 
oblalenoss, which induces a rotation of tho SS orbital 
piano of about 0. 3 dog/hr about tho polar axis. Since tho 
Tug must return to tho SS orbit and the SS does not have 
performance capability to substantially alter Its lino of 
nodes, the Tug is assumed to make all necessary nodal 
corrections for rendezvous. 

As a sample case, payload capability Is dorived as 
a function of the declination of tho outgoing asymptote 
with and without perigee propulsion for a typical cryo- 
genic Tug configuration and an inner planet mission. A 
discussion of optimal trajectory profilos and SS launch 
time constraints Is Included, 

Analysis 

Trajectory Profiles and Assumptions 

A trajectory profile for the perigee propulsion case 
Is Illustrated In sketch (a). Tug and payload are de- 
ployed In a low Earth orbit by tho Space Shuttle (SS). 

The first outbound burn is inltlatod at an optimal point 
along tills orbit. The length of tho bum is determined 
aurlng tho optimization process in solving the two-point 
boundary value problem. At tho end of this bum an 
elliptic orbit Is established with a relatively long period. 
Tug and payload then coast along this orbit to a point 
just before perigee, whore a second outbound burn is 
performed that accelerates the payload to tho given in- 
jection conditions, The timo at which the second out- 
bound burn Btarts is also determined during tho optimi- 
zation. Planetary tnj action conditions are defined by 
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specifying the mi union energy and the declination of the 
outgoing asymptote. The right ascota ion of tile out going 
asymptote is not included as a final condition, because 
ony rlRht ascension ma ■ be obtained by selectinR the 
proper S8 launch time. This Is discussed In more detail 
later. After the Tur reaches the Riven injection condi- 
tions, the Tur thrust is terminated and a fixed coast 
phase is initiated. During this coast the Tur and the 
payload are separated, and the Tur is turned around to 
orient It for the retrobum that follows. DurinR this 
retroburn the enerRy is reduced below Earth- escape en- 
ergy, and the Tug enters an intermediate elliptic return 
orbit. Retrobum terminates when a given total trip 
time is satisfied. Total trip time Is defined as the 
elapsed time from the start of the first outbound bum to 
return to the SS orbit for rendezvous. The retrobum Is 
followed by a coast phase to the apogee of the intermedi- 
ate orbit, where a small perigee correction bum is exe- 
cuted. This is followed by a coast phase to perigee and 
a final rendezvous bum. 

The flight profile for the case without perigee pro- 
pulsion is the sume as Just discussed with the exception 
that the first outbound bum continues until the given In- 
jection conditions are satisfied and thus the coast in 
elliptic orbit and the second outbound bum are elimi- 
nated. 


To olituln solutions to the two-point boundary vulue 
problem and to simplify the anulysls, a number of as- 
sumptions were made. These a»sum|itlon» are discussed 
in the following paragra, .is 

Hi The perigee correction and final clrculurlzatiun 
bums were assumed to lx- ideal impulsive. This as- 
sumption was made to reduce the sensitivity and conse- 
ipienllv Improve the convenience characteristics of the 
two- point boundary value problem. If the corresponding 
state equations are stable, the adjoint equations will have 
unstable roots. Since total trip time for the problem 
under consideration is of the order of I day (based on the 
results of rid. I), errors Introduced In the numerical in- 
tegration of state and adjoint equations will lx- greatly 
amplified. These errors will affect the finite difference 
partial derivatives and consequently the convergence of 
the Newton- Kaphson iteration. With lids assumption, the 
problem is numerically integrated only to the end of the 
retrobum, and the remaining portion of the trajectory Is 
calculated In closed form (Impulsively). 

(2) A circular SS orbit Is assumed. This assumption 
is made to eliminate the constraint on the lint- of apsides. 
For elliptic SS o riil is the Tug's orbit at the end of final 
rendezvous bum, besides being in the same orbital plane, 
must huve Its line of upsides coincident with that of the 
SS orbit. I'sing a circular SS orbit removes this con- 
straint, and rendezvous Is accomplished by small 
changes In total trip time. The trip time has to lx- ad- 
juster! so that the Tug and SS will be at the same point 
along the orbit at final rendezvous I aim completion. 

(it) A spherical nonrotuting Earth model Is used. 

This assumption is made to simplify the equations of mo- 
tion ami the adjoint equations. A nonspherical Earth 
model could be Included with relatively little change In 
the analysis. As a result of this assumption the nodal 
precession of the intermediate Tug orbits becomes zero, 
and the Tug corrects for SS orbit nodal precession only. 
This gives somewhat conservative Tug performance 
at. "o, If the Tug orbits were allowed to process, the 
total nodal correction required of the Tug would be 
slightly reduced, and consequently payload capability 
would Increase. Nodal precession of the Tug orbits for 
the nominal mission is less than 0. 2 deg without perigee 
propulsion and less than 0.4 deg with perigee propulsion, 
as compared with an SS orbit nodal precession of approx- 
imately H deg. 

(4) The SS orbit nodal precession Is competed from 
the following equation: 




This equation was dorlved In reference 2, Nodal pro- 
cession depends on tho BS orbital parameters and Uio 
total trip time T D , 

Bnalc Eqtir'Jons Governing tho Problom 

Varlnblos and otlior notation used In tho following 
dlseuasion aro defined In tho appendix. Kquutlons de- 
scribing Uio flight of a rocket In an Invorso-squaro 
gravitational Hold nro 


G - + V#> f 

r 3 m 

(2) 

?= V 

(3) 

m » -/!(£) 

(4) 


<G) 

f.f- l«o 

(0) 


In those equations, r, V, and m represent Instan- 
taneous radius, velocity, and mass, respectively; V Q 
Is tho engine oxhauBt velocity; /? is tho mass flow rate; 
and f is tho unit thrust direction, which must bo deter- 
mined, Tho state variable £ was Introduced to romovo 
tho explicit llmo dependence of mass flow rato from Uio 
equations of motion so that the Maximum Prlnolplo could 
bo applied to this problom. Tho superscripts - and * 
represent voctor and unit vector quantities, and • Is 
tho total derivative of tho particular varinblo with re- 
spect to ttmo. 

By using tho Hamiltonian formulation of tho varia- 
tional calculus, tho costalo equations may be obtained 
from the so-called Hamiltonian on each subarc. The 
Hamiltonian Is 

Il=l-V+p-r +o-m+r£ +y(f-f- 1) (7) 

In this equation, X, J7, o, and t aro tho adjoint varia- 
bles associated with tho problem. Tho constraint asso- 
ciated with tho thrust direction Is also adjoined to H by 
y. From equation (7) the time derivatives of the adjoint 
variables aro given by 



(8) 

? = — [X- 3{X • r)r] 
r 3 

(9) 

V„/3 

s-tx- 0 

(10) 

m 2 


r=-K^ 

(U) 


0 £ 


Ko-J2(X.Q-cr (12) 

m 

Tho thrust direction that minimizes tho Hamiltonian Is 
given by 

f«X (13) 

and 

V.tf 

V> \ (14) 

2m 

State equations (2) to (5) and adjoint equations (8) to 
(11) make up a sot of first-order nonlinear differential 
equations that must bo numerically integrated In order to 
obtain a solution to tho two-point boundary value problem. 
Optimum thrust direction Is along Uio unit vector given 
by equation (13). 

Trajectory Constraints 

In tills section, constraints on 'ho problom are dis- 
cussed In more detail. Sketch (b) Illustrates a simpli- 
fied trajectory profile from a given Initial orbit to the 
ond of rotrobum. Events a ( and b ( Indicated on the 
sketch are Uie times at tho beginning and ond of phase 1, 
respectively, 



The first outbound bum Is Initiated at time a 1 
(Oj — 0) at an optimal point In llio initial SS circular In- 
clined orbit. Main engine thrust terminates at time bp 
which is selected to give maximum payload. At this 
time on elliptic orbit is established and the Tug coasts 
along this orbit to a point near perigee (time b 2 ). Sec- 
ond outbound bum starts at time a 3 (a 3 = b 2 ) and con- 
tinues until the specified payload Injection conditions are 
reached at time b 3 . Time a^ (a^ = b 3 ) Is tho beginning 


where 


(15. 25) 


oC n fixed payload- separation coast that continues until 
time b 4 . Tug retroburn starts at tlmo a s (u c « b^) and 
ends at bg. Rotroburn Is followed by Ideal impulslvo 
npogoi and porigoo buns and associated coast phases. 
With lids tlmo convention, constraints q, governing the 
problem can bo written mathematically as follows! 


qj a via,) - r x = 0 

(15. 1) 

q 2 = V(n x ) - = 0 

(15.2) 

q 3 = r(ai).V(ai) = o 

(15. 3) 

q,j = h(nj) . t 3 - cos I 4 = 0 

(15.4) 

l B = -hjttj) • f 2 - cos Rj sin ij = o 

(15.0) 

q 0 = mfa^) - = 0 

(15.0) 

“£(«!> “ a i " 0 

(15.7) 

q g = r(a 2 ) - r (b x ) = 0 

(15. 8) 

qg = V(02) - V(b x ) =0 

(15.0) 

(t 10 = m(a 2 ) - m(b 4 ) = 0 

(15.10) 

q^l = £ (a 2 ) “ £ (b 4 ) = 0 

(15.11) 

“l2 a ^ a 3 } " ~< b 2 ) a ° 

(16. 12) 

qi3 = V(a 3 )- V(b 2 ) = 0 

(10. 13) 

“14 “ ™{ r 3) “ m &2) ~ 0 

(16. 14) 

fl :5 = t(a g ) - £ (b 2 ) = 0 

(15, 15) 

“io “ *< a 4> " *< b 3> n ° 

(15. 10) 

“17 = V ( a 4> " V(b 3> “ ° 

(16, 17) 

q lg = m(n 4 ) - m(b 3 ) = 0 

(15. 18) 

qi 0 “£ta 4 ) - £(b 3 ) =0 

(15, 19) 

q 20 = E[r(b 4 ), v (b 4 )]- E D =° 

(15. 20) 

q 2 i = V>[r(h 4 ). V(b 4 )] - v> D = 0 

(15.21) 

q 22 = r(a 5 ) -r(b 4 ) = 0 

(15. 22) 

q 23 = V(a 6 ) - V(b 4 ) = 0 

(15.23) 

q 24 = ^ ^ (^ 4 ) “ ® 

(15.21) 



q 25 «r(b 5 )-l) o < r D > D0 


“27 ' 


- V(b 0 ) • * 0 

(15.20) 


(15. 27) 

“28 R n 2 - b l “ 0 

(15,28) 

“20 “ a 3 _ b 2 “ 0 

(15.29) 

“30 D a 4 " b 3 “ 0 

(1G. 30) 

“31 n a G “ b 4 ” 0 

(16.31) 

2 = a 4 + At - b 4 = 0 

(1G, 32) 

£(a 4 ) +At- £(b 4 ) =0 

(15, 33) 


Equations (15, 1) to (10.7) express fixed initial conditions 
at Tug departure from the Initial SS orbit. Equations 
(15. 8) to (15. 10) givo tho continuity of radius, velocity, 
mass, and tho dummy variable £ at tho end of tho first 
and second outbound bums and at the beginning of tho 
second outbound bum. Equations (ID. 20) and (1G. 21) 
oxprcsB tho constraints associated with a planetary mis- 
sion, namely, at InjocUon tho payload must have a given 
energy and declination of tho outgoing asymptote. Equa- 
tions (1G, 22), (15.23), and (1G, 2'1) give tho continuity of 
position, velocity, and £ at injection. Satisfying equa- 
tions (15,26) and (1G,2C) assuros that tho Tug Is In an 
orbit Ihnt has tho snmo Inclination at rctrobum comple- 
tion qb tho SS will have at rendezvous. In other words 
the flight Is planar following tho retroburn, Tho vector 
h Q in equations (16. 25) and (15. 20) is tho unit angular 
momentum vector of tho precossed SS orbit to which tho 
Tug must return, Tho vector h Q is a function of tho 
total trip ttmo. Equation (15, 27) is tho constraint on 
Tug mass required at b g to perform tho apogco and 
perigoo velocity changes and return the Tug empty to 
tho final SS orbit. Equations (15.28) to (15.31) express 
the continuity of tlmo between different trajectory 
phases, equations (15.32) and (16, 33) are tho fixed 
payload- separation coast time constraint and tho con- 
straint on £, respectively. 

With those constraints the augmented Hamiltonian 
(H + ) to bo minimized can bo written as follows; 

5 33 

H* = m(a 6 ) - m(b 4 ) + £ Hj + Y, e i“i {1B) 
1=1 1=1 


where m(b ( j) - m(n g ) represents tho payload weight 
dropped during tho coast phase, which is to bo mrxi- 
mized; and IIj nro Hamiltonians during tho two out- 
bound bums, the two coast phases, and the retroburn. 



Constraints q^ oro adjoined to H* by using arbitrary 
multipliers cj. ouxlliary vnrlntlonnl boundary co> '>i- 
tlons may now bo dorlvod from tills augmented Hamll- 
tontan. 

Auxiliary Vnrlatlonnl boundary Conditions 

Tho following auxiliary variational equations must 
bo satisfied by using tho Maximum Principle and tho 
augmented Hamiltonian In equation (10) In order to max- 
Imlzo payload; 

X(<V ° + ^(“l) + c -l v v fl -l + f D V V q G (17< *> 

p(a 1 ) « e 1 r(n 1 ) + c,V(a t ) + e 4 V-q 4 + c c V-q s (17.2) 


m(b,.) 

, 0( b 5 )_- r ^v ? AV T 

* n 


*JSL\ 

r(b G ) 1 

[‘V r(b rj )]r(l J G )} “0 

(17,20) 


lHftj) - IHty 

(17.21) 


H(n b ' « U(b 2 ) 

(17.22) 


X(n 2 ) "X(b 4 ) 

?(a 2 ) “f»(bj) 
o(n 2 ) = c(bj) 

T( a 2 ) 

X{Og) “X(b 2 ) 

p(a 3 ) - f»(b 2 ) 
o(a 3 ) = a(b 2 > 

T(a 3 ) °T(b 2 ) 

M® 4 ) =X(b 3 ) 
q (a 4 ) = p(b 3 ) 
tf(a 4 ) c cr(b 3 ) 

X(a 6 ) =X(b 4 ) + c 20 V^E(b 4 ) + e 21 V^(b 4 ) 
p(a 6 ) =?(b 4 ) + e 20 V r E(b 4 ) + £ 21 V r v>(b 4 ) 
cf(a 5 ) = 1 
o(b 4 ) = 1 

T (a 4 ) - y(b 3 ) +r(a 5 ) - T (b 4 ) = 0 
X(b 5 ) + [c 25 r(b r) ) + e 20 V(b 5 )] • VyT 
m(b-) 

+ ‘ T(b 5) ~~~ V V AV T 
v o 


(17.3) 
<17. 4) 
(17. B) 
(17.0) 
(17. 7) 
(17. 8) 
(17.9) 

(17.10) 

(17.11) 

(17. 12) 

(17.13) 

(17. 14) 

(17.15) 

(17. IB) 

(17.17) 

(17. 18) 


(17. 19) 


II(a G ) + H(a 4 ) - ll(b 3 ) - ll(b 4 ) * 0 (17. 23) 

»( b G> “ [e 2 G* r < b 5) + C 2Q V(b c )] • ^ (17. 24) 

r(b 0 ) = 0 (17, 25) 

Equations (17. 1) and (17, 2) contain six equations and five 
arbitrary constants. Therefore, there is a variational 
condition that tho Initial multipliers must satisfy . ■ max- 
Imlzo payload. Equations (17, 3) to (17. 13) gtvo tho con- 
tinuity of adjoint variables at tho end of tho first and 
second outbound bums and at tho beginning of tho second 
outbound burn. These equations eon bo satisfied directly 
by setting tho multipliers at tho beginning of those phases 
equal to the corresponding multipliers at tho end of the 
previous phase. Multipliers at the end of tho payload- 
separation coast aro discontinuous, as shown by Dera- 
tion,? (17, 14) and (17. ID). This discontinuity w! V. 
computed later. Multiplier a Is continuous an < ■ ■ 
to 1 at tho end of tho payload- separation coast pin v 
(oqs. (17. 10) and (17.17)). Equations (17, ID) and (17.20) 
must be satisfied at retroburn completion, Thoy contain 
six equations and three unknowns (c 2g , e 2( ,, and o(b g )), 
Tho three variational conditions available for optimiza- 
tion will bo evaluated lator. Tho Hamiltonian Is contin- 
uous at the completion of the first outbound burn and tho 
start of tho second outbound bum, as given by equations 
(17. 21) and (17. 22). Equation (17. 23) gives the relation 
the Hamiltonian must satisfy at tl>o end of tho second 
outbound bum and tho start of tho retroburn. Gradients 
of return time T, total delta velocity AV.j, - AV a * 

AV , and declination ip wore derived in reference 1 
anti will not be reproduced horoln. Equation (17. 24) 
gives tho Homiltonton at retroburn completion, Tho 
vector h Q Is a function of tho total trip time, or return 
time, and the required partial derivatives will bo evalu- 
ated lator. Return time is defined as tho time from 
rotrobum completion to return to tho SS orbit. Equa- 
tion (17.25) is tho boundary condition on tho multiplier 
associated with tho dummy variable £ . 
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Variational Finn! Condition at Departure 
from Initial SS Orbit 




( 20 ) 


V 




To dorlvo tho variational final condition contained In 
cquatlona (17. 1) and (17. 2), dcflno tlio vector 

C(t> = X(t) x V(t) + j?(t) x r(t) (1 8) 

whoro x represents tho vector or croas product In this 
equation, Tho vector C if, a constant of motion. This 
can bo shown by tnklnq tho timo dorlvattvo of C and 
substituting equations (2), (3), (8), and (0) for ' T, 

and Ui roapoetlvoly, Tho ttmo derivative ut ( i , 

out to bo Identically zero; thoroforo, C is a co. 
motion. 

Tho variational condition at doparturo from tho Ini- 
tial orbit may now bo computed as follows. Compute tho 
C vector at a^ giving 

From this equation tho variational final condition 1 b ob- 
tained by taking tho scalar product with lifa^), giving 

0(3^ • itOtj) = 0 (20) 


x V(b 4 ) + V- r v> xf(b 5 )] • fi Q 

Given Cgi an(1 X(b 4 ), lot 

X“X(b 4 ) + e 2l V?V V' i 

Will) this definition, equation (17, 14) becomes 

X(n c ) « c 20 V(b ( j) + X (27) 

To determine tho discontinuity tn X and /7, c 2Q mu °t 
bo evaluated. This may be done as follows. From equa- 
tion (11) 

T (a fl ) - T (b,j) - -d(a e )K(a 5 ) (28a) 

and 

r(tt 4 ) - r(b 3 ) - /?(b 3 )K(b 3 ) (28b) 

From equations (28) and (17. 18) the following equation is 
obtained; 

/J(a c )K(a c ) = /3(b 3 )K(b 3 ) (20) 


Evaluation of Discontinuity In State Variables 
at End of Payload-Separation Coast Phaso 

To ovaluato the discontinuity In X and JT at timo 
a s , two arbitrary constants (c 20 and e 21 ) must bo eval- 
uated in aquations (17. 14) and (17. ID), Those constants 
aro caluatod In this section. Take the vector produots 
of equation (17. 14) with V nnd (17, 15) with T, add the 
resulting vector equations, nnd apply equation (18), 

This proooduro will result In 


Since a => 0 on a coast phaso (ff= constant), equations 
(17.13), (17.10), and (17. 17) gtvo or(a g ) = e(b ;j ). With- 
out loss of generality, assume )3(a g ) o /j(b 3 ) 4 0, Then 
equation (20) reduces to 

m(a r ) 

< 3 °) 

m(b 3 ) 

Now X(a g ) may bo substituted from equation (27) and tho 
resulting equation solved for e 20 , giving 


C(a g ) - C(b 4 ) = c 21 [Vy<p x V(b 4 ) + V-ip x r(b 4 )] (21) 

Since C Is constant during each phaso 

C(a G ) = C(b g ) (22) 

Using equations (17.19) and (17.20) yields 

°(fc 5 ) = \ X [«26^ + *25* r < b 5>] < 23 > 

Tilts equation was derived by using tlio results of refer- 
ence 3, whoro it ts shown that Vyyx V + V^y x r = 0 
if y = y(r, V, r-V). Since the return time T and tho 
apogee and perl jeo doltu velocities aro of tills form, tho 
corresponding vector product sums are zero. From 
equation (23) 

C(b g ).h o =0 (24) 

Taking tlio scalar product of equation (21) with h Q Uien 
gives 


e 20 


. .1 

V(b 4 ) 


X-V(b 4 )±l/[x.V(b 4 )] 5 


m(a g )Mb 3 ) 


i*l 2 


L m < b 3> J 



(31) 

Choice of + or - sign can bo determined as follows. 
Assume the payload Is zero, that is, m(a g ) = m(b 3 ). 

Also If the declination Is not specified, equation (20) be- 
comes X = \(b 4 ), Furthermore, If payload- separation 
coast phase Is zero, X(b 4 ) = X(b g ), With these assump- 
tions, equation (31) becomes 

*2° “ ' “ * V< b 4) ± *< b 4> ‘ V( b 4>] 

A choice of - sign Implies that e 23 a 0 and that the 


G 



thrust direction is continuous ill a B und Is essentially 
pointed along tho velocity vector. Howovor, to bring die 
Tug back, velocity must bo reduced bolow osenpo veloc- 
ity and tho thrust direction should bo opposite to tho ve- 
locity vector, Therefore, tho + sign should bo used In 
equation (31), 

Variation!!! Final Conditions at 
Itotroburn Completion 

In this section, variational final conditions will be 
dorlved at retorbum completion. Substituting equations 
(17. 10), (17,20), and (17. 25) Into equation (17. 24) gives 


f 20 : 


X(b 0 ) ■ h(b f) ) 


V(b c ) 


(31) 


and 


C 2G 


H(bfi) • filb„) 
fi(b c )< n Q 


r(b Q ) 


(30) 


Two variational conditions may now ho obtained by taking 
tho scalar product of equation (17, 10) with Vy AV,j. and 
equation (17, 20) with v-^ AV T and using equation (33): 


- [c 25 r(b 0 ) + e 20 V(b B )] • 



^X{b c ) +X(b 5 )Vy AV, r + ““ {>V[V V(b c )]v(b c )) 


x 


— £- V V T . r(b B ) + V-T . V(b 5 ) 

r 3 < b G > 


- u(b E ) 


m(b s ) 

v 0 



r 3 (b fi ) 


VyAV T > r(b c ) +V-AV t - V(b B ) 




V< b D> 


L m(b c ) 


- u(b B ) 


0b„ 


“ [ e 25 r(b 5 ) * «20^< b 6)] ’ TZ 


(32) 


(30) 


+ MV y r AV T ft, - Po ' ?(b s )]r( b 6 )) 

+ [wW + <2G<M ~ V 1 ) * V r AV T = 0 


(37) 


However, 


In these equations, afi Q /0T D is given by the vector 


- -2- VyT. r(b 6 ) H-v-T. V(b 5 ) 0^2 
r 3 (b 5 ) (lt 



(38) 


and 


whoro 


- ^ A V + V r AV T • V <V = ~ 

r 3 (b c ) dt 


jjjo 

ot d 


_ ea 


&T 


D 


(r n v i3 


r 13 V ll’ r 12 V 13 " r 13 V 12’ 0) 


These equations are valid on a coast phase, whoro tlio 
delta velocity required to return the Tug to the SS orbit 
remains constant; and d AV,j,/dt = 0. Similarly, tho 
change in return time Is equal to tho change in time 
nlong tho coast are but opposite In sign, that is, 
dT/dt = -l. Substituting these results Into equation (32) 
givoe 


u(b 5 > 


V< b G> 

m(bg) 


(33) 


Tho components of r^ and in tills equation are 
givon by 


r lj = r l l£ ] 


Vy-V'j 


Tlio tldrd variational condition is given by using equa- 
tions (17.24), (17. 2G), and (33) as 


The arbitrary constants e 2B and ego can bo evaluated 
by taking tlio scalar product of equation (23) with r(b B ) 
and V(bjj), giving 
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— 2— x(b ( ) • r(b 0 ) +P(b 5 ). V(b 0 ) 

r 3 < b 0> 


„L 


- [<21>G> + f 2fl^< b D>] 1 " 0 < a0 > 


Variational ConilitlonH nt Completion of First Outbound 
Rum anil Start of Second Outbound Burn 


Integrating equation (11) and applying equations 
(17.0) and (17,10) gives 


anti 




m(b x ) 


o(a 3 ) 


,w 

in(n ;1 ) 


(<10Q) 


(-10b) 


Since a is constant on a coast phase, equations (17, S) 
ami (17.9) givo o(a 3 ) = o(b 4 ). Therefore, slnco mfbj) = 
m(i j), one vaiintionnl final condition is given by 

X(a a ) - X(b x ) e 0 (-11) 


Tho second variational condition is derived based on the 
constancy of tho Hamiltonian on each phase and equntlons 
(17,22) to (17,25). It is gH . n by 

- — S — X{lf 2 ) < r(b 2 ) +Ji{b 2 ) • V(b 2 ) + T (bg) 
r 3 (b 2 ) 




(42) 


whoro 


r(b 2 ) = 



K^dE 

DC 


+ r* Kifldt 

L * 


(43) 


If tho Tug has a constant flow rate, r(bg) is equal to zoro 
and r need not bo numerically integrated, 

Initial and Final Conditions 


Initial 

Final 


conduit, ( 

conditions 


if (u t ) 

‘r(b c )-f. o -=0 

(44.1) 

V r (a x ) 

V(b 0 )’ b o 0 

(44.2) 

6(aj) 

tf[r(b,,), V(bj)j- - 0 

(44,3) 

6(Uj) 

liquation (30) 

(44.4) 

u(a x ) 

Equation (37) 

(44, 5) 

m(n c ) - m(b|) 

*V V D 

m(b fl ) - nijjO - 0 

(44.0) 

t D 

Equation (30) 

(44.7) 

b l 

Equation (4l) 

(44. 8) 

b 2 

Equation (42) 

(44. 9) 

Mn x ) 

o(a 5 ) *= a(b 4 ) = 1 

(44.10) 

Xfaj) 

C(a x ) * hfaj) = 0 

(44,11) 

b 3 

E[r(b 4 ),V(b t )]- E d -0 

(44, 12) 

b 5 

b G -a 1+ T-T D =0 

(44. 13) 


Final condition (44,10) can bo satisfied by scaling 
tho problem by Ma x ). Thoroforo, It is omitted from the 
numerical iteration. Final condition (44,11) is used to 
oompulo ^(a x ), as shown in reference 4, Final condi- 
tions (44, 12) and (44, 13) aro used as cutoff conditions to 
determine b 3 and bg, With those simplifications the 
number of final conditions lo be satisfied for the perigee 
propulsion enso is reduced to nine, namely, conditions 
(44. 1) to (44. 9). For tho case without perigee propul- 
sion, conditions (44,8) and (44,9) do not apply, and tho 
number of final conditions to be satisfied Is reduced to 
sovon, 

Tho two-point boundary value problem Is solved by 
using a Nowton-Raphson numerical iteration technique, 
Tho iteration is terminated when tho percent of predicted 
payload change and the absolute value of the normalized 
final conditions aro less than 10“ e . 


Tho inlt* . adjoint variables aro unknow and must 
be guessed at in solving the two-point boundary value 
problem. As was shown in reforenco 4, multipliers X 
and p can bo computed from physically more meaning- 
ful parameters, such as the vehicle pitch attitude i p, 
pitch rate ip, yaw attitude 5, yaw rate &, and X(a x ) 
and ^(a x ). With those changes in variables the initial 
and associated final conditions aro given by 


Results and Discussion 

Tho method and equations derived in tho previous 
sections were appliod to n typical cryogenic Tug config- 
uration and a sot of planetary Injection conditions. The 
rosults of this analysis aro presented in tills section. 

Propulsion and weight characteristics for the re- 
usable Tug aro given In table I. The initial wolght of tho 
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Tug and payload Is assumed to t )0 fixed to remain ’.vlthtn 
SS paylond capability. To satisfy tills constraint, Tin; 
propellant load must bo varied with payload weight to 
maintain tills Initial weight, Engine thrust corresponds 
to tho nLr-!0 engine currently In use on the Centaur 
cryogenic uppor stage, but engino Bpeclflo impulse re|>- 
rosentB a slnto-of-lhe-nrt uprated version of tho same 
engine. To oocoiuit for losses and flight performance 
reserves, ongbic specific impulso is reduced by 2 per- 
cent, 

Tho sample case selected (table il) Is a typical Inner 
plnnot (Mnrs, Venus) mission, with a mission energy of 
12 (km/seo) 2 , an Initial SS elreular-orbll altitude of 
185 km, and n SS orbital Inclination of 28. 5°, Tills SS 
orbital Inclination was chosen based on tho results of 
refcrenco 1, whoro It Is shown that maximum Tug per- 
formance Is obtained (for declinations less than 30°) 
when departing from this SS orbital Inclination, 

Payload- separation coast time (tho time from outbound 
burn cutoff to tho start of retrobum) Is assumed to bo 
10 minutes. 

To show the payload gain duo to perigee propulsion, 
Tug payload capability without perigee propulsion Cor op- 
timum trip times Is also presented. 

Figures 1 lo 7 present rosults without perl geo pro- 
pulsion; figuros 8 to 10 contain results for porlgco pro- 
pulsion oasos, Thfr data aro given as a function of tho 
declination of tho outgoing asymptote (DLA), which was 
parametrically varied. 

It Is also shown In reference 1 that Ihoro aro two 
solutions for oach declination and that tho paylond curves 
aro symmetric about a 0° declination, Thereforo, tho 
data presented herein nrc restricted to posltlvo deollna- 
ttons. night ascension of tho outgoing asymptote was not 
specified In tho annlysls, and consequently an "optimal 
psoudo right ascension" Is generated for each DLA, 
Pseudo right ascension is dcflnod as tho longitude of Jio 
projection of tho outgoing asymptote in tho equatorial 
piano measured from the Initial SS orbital ascending 
node. Right ascension may bo computed by adding the 
psoudo right asconslon to Uio longltudo of tho SS orbital 
ascending node measured from tho vernal equinox. 

Since tho longltudo of tho SS orbital ascending node Is 
determined by SS launch time, any desired right ascen- 
sion may bo obtained by selecting tho proper launch 
time. 

Optimal pseudo right ascension 1 b presented In fig- 
ure 1 for die case without porlgco propulsion. To dis- 
tinguish between Ihe two solutions on all figure*, a solid 
and a dashed lino aro used, Tho solution given by tho 
dashed lino will bo referred to as solution 1, and tho ono 
given by tho solid lino as solution 2. In tho region of 
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single solutions, where DLA exceeds SS orbital Inclina- 
tion, a solid lino Is used. Note Unit for 0° DLA, tho out- 
going n'V .’iplotos associated with the two solutions point 
InoppoeAo directions, As DLA in Increased the outgoing 
asymptotes movo closed together until the two solutions 
degenerate Into a single solution, At this point the out- 
going asymptote has a pseudo right ascension of approxi- 
mately 00°, For DLA' s beyond tho SS orbital inclina- 
tion, tho pseudo right ascension remains nearly constant 

Tho Tug departure point frev tho Initial SS orbit was 
determined during tho optimization, and the results arc 
presented in figure 2, Departure arguments of latitude 
of tho two solutions aro also 180° apart for 0° DLA and 
converge to about 250° as Uio DLA approaches tho SS 
orbital Inclination, 

Optimum total trip lima (tho llmo from Tug depar- 
ture from tho Initial SS orbit to return to tho rendezvous 
orbit) is given In figuro 3, Total trip times for the two 
solutions aro nearly tho samo and decrease from approx- 
imately 20 to 21. 0 hr us DLA Is Increased from 0° to 
30°, Total trip time vanes Inversely with the duration 
of Uio retrobum, Tho shorter tho retrobum, the more 
elliptic tho return orbit, and consequently Uio total trip 
time increases. Velocity loss during retrobum is pro- 
portional to Uio burn time. Therefore, as total trip time 
Is Increased, tho velocity loss during retrobum de- 
creases. llowevor, the nodaL correction duo to SS or- 
bital procoBSion will Increase with total trip time, which 
increases tho velocity loss to correct for ndal preces- 
sion. Tliorcforo, tlioro Is a balance between those two 
effects that determines tho optimal trip time. Variation 
in trip ttmo Is duo to changes In oruttal geometry and 
associated nodal correction during outbound and retro- 
bums, Tho sharp Increase in total trip time for DLA's 
larger Uian the SS orbital Inclination can bo explained 
when changes In orbital Inclination during outbound bum 
aro considered, it Is woll known from orbital geometry 
that orbital Inclination nt Injection must bo at least as 
lnrgo as the magnitude of tho desired DLA, Thereforo, 
largo Inclination changoB nro made during boUi outbound 
bum and retrobum for DLA 1 s Inrgor than Uio SS orbital 
Inclination. Departure arguments of lnlitudo for these 
casos aro approximately 200°, and Uio outbound bum 
covers about a DO 0 aro. Therefore, Uio complete bum 
occurs before nodal passage. Since orbital Inclination 
must be increased, ns was mentioned earlier and the 
lino of nodes tends to move retrograde during Uio In- 
clination change, it becomes easy to make largo nodal 
changes during tho outbound burn. Those nodal changes 
aro in tho some direction as tlio SS nodal precession, so 
longer trip times bocome optimum. 

Nodal correction during Uio outbound bum and total 
nodal correction (during outbound bum and retrobum) 
aro given In figure 4. Total nodal corrections for tho 



two flotations tiro vory close and differ only beenuso of 
tha slightly different trip times for tlio two solutions. 
However, the outbound nodal correction for one solution 
Is largo and for Iho other It In small, This is explained 
by tlio location of Iho outbound burn. For the case with 
largo nodal correction during outbound burn, tlio out- 
bound bum starts near or after nodal passage (fig. 2). 
SIneo Iho outbound bum arc is approximately 00°, tho 
Tug bum coalinues through maximum latitude, which Ib 
an optimal region In which to change tho lino of nodes, 
For Iho caso with small nodal correction during outbound 
bum, tho hum genornlly tnhes place near tlio ascending 
nodal passage, where it Is difficult to chango tlio node. 
Also noto that for largo DLA's tho total nodal correction 
Increase!! faster than the nodal chnnga during tho out- 
bound bum; that Is, larger nodal changes are mado dur- 
ing tho rotrobuni, This Is expected, since largo Inclina- 
tion changes arc being nindo during HiIb bum to return 
tlio Tug to tho SS, and it becomos optimum to comblno 
nodal correction with tho Inclination change. 

Figure 0 presents tho inclination of iho hyperbolic 
orbit containing tlio outgoing asymplolo. Thoro is a 
maximum inclination chango during outbound bum of 
41, D dog for DLA's loss than the SS orbital Inclination, 
For tlio first solution (dashed lino) tho outbound burn 
starts ahead of tho ascending lino of nodes (fig. 2) and 
ends post It, In tills regton ll 1 b difficult to change tho 
lino of nocJos, and nodal change (fig. 4) is Bmall, To 
movo the node retrograde for those solutions, tho or- 
bital Inclination must bo increased, ns oliown In figure G. 
For tho second solution (solid Uno) tho outbound bum 
starts past nodal passage (fig. 2) and continues through 
maximum latltudo, which Is an optimal region la which 
to chango tlio nodo. Most of tho nodal chango Is done 
during this bum (fig, 4). Slnco tho bum is past tho node, 
orbital Inclination roust bo decreased to movo tho node 
retrogrado (fig. G). As tho DLA approaches Iho SS or- 
bital Inclination, the Inclination at Injection must bo In- 
creased to moot tlio DLA requirements, and for tlio sin- 
gle solutions tho orbital Inclination Is Increased to coin- 
cide with tho DLA, 

Tug payload capability without perigoo propulsion 
far Iho two solutions is given in flguro 0. Tho two solu- 
tions hnvo nearly tho same porformnneo (maximum dif- 
ference, 20 kg out of 47G0 kg). Payload capability de- 
creases with Increasing DLA. This decrease Is ap- 
proximately 00 kg between 0° and 30° DLA, For DLA's 
higher than 30° tlio payload decreases rapidly because 
of tlio largo Inclination changes required, as discussed 
oarlior. 

Velocity losses encountered durin b outbound and In- 
bound trajectory portions ara presented In figure 7, 
Voloelty loss Is defined as tho difference between tlio 
velocity supplied by the Tug (based on tlio amount of 


propellant used and the Ideal rocket equation) and the 
ideal mission Impulsive velocity, Ideal mission Impul- 
sive velocity Is defined as tlio velocity needed above 
circular-orbit velocity at n 183-km altitude In order U> 
reach a given energy, Total velocity loss gradually In- 
creases with Increasing DIJV for liiA's less than .30°. 

As DLA la Increased beyond 30°, the velocity losses in- 
crcaso very rapidly, again because of tlio large Inclina- 
tion changes being made, Noto Dial the solutions given 
by the dashed Uno have larger total velocity losses than 
tho solutions given by the solid line, yet tlio payload for 
these solutions Is higher (fig, 0), This apparent Incon- 
gruity Is duo to tlio unequal effect of velocity loss on pay- 
load during outbound and Inbound trajectory portions. 

Tho effect of tho outbound velocity lass on payload Is al- 
most twice the effect of tlio Inbound velocity Iosb, 
Therefore, for tlio solutions given by tho dashed lino the 
payload Increase duo lo lower oullxiund velocity losses is 
greater than tlio payload decrease duo to higher Inbound 
velocity losses when compared with tho solutions given 
by tlio solid lino, 

Tho outbound velocity losses of figure 7 Indicate that 
large velocity losses nro associated with tho ftnlto thrust 
level of tlio reusable Tug. These losses can he lcdueed 
substantially by using Die tcclmiquo called "perigee pro- 
pulsion, " With porlgco prape'slon the outliound bum Is 
split into two separate bums. The first outbound bum 
puts tho Tug and payload Into an elliptic orbit, Tlio Tug 
Is allowed to coast almost to the porigeo of this orbit, 
where a occond outbound bum takes place. The Bccond 
outbound bum continues until tho target conditions aro 
satisfied. After payload Injection nt tlio specified condi- 
tions, tho flight sequence Is tho same as that without 
perigee propulsion, Tho lengths of the first outbound 
bum and tlio coast phoso following It arc optimized, dial 
Is, they nro selected lo maximize payload. 

In figures 8 and 0, the pseudo right ascension and 
the departure argument of latltudo are given, respec- 
tively, as a function of DLA for the caso with perigee 
propulsion. These, figures aro vory similar to thoao 
discussed for tho ease without porlgoo propulsion and 
will not bo elaborated on further. 

Optimum lolat trip time Is presented In flguro 10. 
For those eases the first solution has a shorter trip 
time (dashed lino) than tlio second solution. The first 
solutions also havo shorter first outbound bums, as 
shown in figure 11, Tlio first outbound bum varies with 
varying DLA from about 030 sec to 740 see, as com- 
pared with tho total outbound bum of approximately 
1100 f.ee for die non- perigee- propulsion ease. Tho pe- 
riod of die elliptic orbit achieved nt first outbound bum 
completion varies between 4 and 0 hr, Tho second out- 
bound bum is Initiated along Uiis orbit at die optimum 
true nnonmllcs given In flguro 12, 
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Tito total outbound nodal correction am) the Inclina- 
tion of the hyperbolic orbit containing the outgoing nu- 
yinptoto nro presented In figures l.'l and 11, resjiecUvely, 
and nro similar to those without perl poo propulsion. The 
larger deference In total required nodal correction be- 
tween the two solutions with porigeo propulsion than be- 
tween the two solutlona without perigee propulsion iB 
caused by a larger dlffereneo in total trip lime. 

Payload capability with porigeo pi’opulnton lu given 
In figure 10. Tho maximum payload difference between 
tiio two solutions is about 100 kg out of about G200 kg. 
Payload agnln drops off rapidly ns DLA Is Increased be- 
yond 30°, Velocity losses with porigeo propulsion arc 
substantially lower than without porigeo propulsion, nu 
soon by comparing figure 10 with Hguro 7, This results 
in tho Improved porformnneo using porigeo propulsion 
neon In figure 1G, 

For comparison, tho Ideal Impulsive payload capa- 
bility without nodal correction and tho payload capability 
including nodal correction with and without porigeu pro- 
pulsion nro given In figuro 17, Ideal payload capability 
for tho baseline Tug analyzed Is G700 kg, as compared 
with a maximum payload capability with porigeo propul- 
sion of G100 kg mid without porigeo propulsion of 1700 kg. 
Tills allows that by using porigeo propulsion tho payload 
loss duo to gravitational and othor trajectory shaping ef- 
fects can be reduced by about 10 percent, Of course, 
doing so requires on additional bum with Its associated 
startup and shutdown losses, ns well as additional guid- 
ance requirements that must bo considered In evaluating 
tho advantages of porigeo propulsion, 

Summary and ConoliiBlons 

Equations aro derived In Hits paper that can bo used 
to maximize tlio payload capability of a reusable Tug at 
energies beyond Earth csenpo With and without tho uso of 
porigeo propulBlon, The analysis includes correction 
for SS orbit nodal precession, while tho total trip time 
(time from Tug departure to return to tho SS orbit) Is 
optimized, 

Based on tho results presented for the bnsoltno Tug 
and baseline mission, tho following conclusions con bo 
made: 

1, For trajectories not using porigeo propulsion tho 
payload lossos duo to gravitational and othor trajectory 
Bhaplng constraints arc large (070 kg out or an Ideal pay- 
load capability of G7G0 kg), 

2, Tho payload loss can be reduced by approximate- 
ly 10 percent by using perigee propulsion, (Payload 
loss with porigeo propulsion Is G70 kg,) 

3, Optimal total trip time varlCB from about 21 to 


20 hr without perigee propulsion and from 1H to 21 hr 
wllli perigee propulsion for a declination range of 0° to 
30°. 

1, Any declination between -30° and *30° may be 
reached by departing from a 2h, 0° Inclined SS orbit, Tho 
variation In payload for this declination range 1» small 
compared with the nominal payload capability, For high- 
er declination!) the SS orbital Inclination muBt bo In- 
creased to avoid largo payload degradation, 

G, Optimal round-trip reusable Tug trajectories in- 
cluding tho S3 orbit nodal precession to energies beyond 
Earth escape with and without perigee propulsion can be 
obtained by a straightforward Newton- Itnphson Iteration 
techniquo uBlngthc methodology developed heroin. 

Appendix - Symbols 

a ttmo at beginning of a phase, sec 

b time at end of a phase, sec 

C defined by eq, (16), kg-Sce 

DLA declination of outgoing asymptote, deg 

E energy, nT/oec 2 

□ eccentricity 

f thrust direction deHncd by cq. (13) 

G gravitational constant of Earth, m'Vscc 2 

H Hamiltonian, kg 

II* augmented Hamiltonian, kg 

h angular momentum, m 2 /sce 

1 Inclination, rad 

O 

J oblatencss parameter (1,021x10 ) 

K kappa function defined by eq. (12), sec 

£ { g rlght-hnndad CurtcBlan coordinate system 

(vector £ j points to tiio initial SS orbital 
Uno Of nodes and £ n points to tho North 
Polo) 

m mass, kg 

P orbital period, sec 

p somilatus rectum, m 

qj constraints 

Il E radius of Earth (G. 37810x10° m) 

r radius, m 

T return time, sec 

Tp total trip lime, see 

t time, sec 

At payload-separation coast time, sec 


, f ,?'AG51 £3 
} ' F £OOB QUMsrfl 
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argument of latitude, rad 

velocity, m/soc 

ai»geo delta velocity, iit/scc 

perigee delta velocity, m/sce 

total delta velocity (AV T • AV a + AV )t ), m/see 

defined by eq, (20), kg-uec 2 /m 

dummy variable 

mass flow rato, kg/sec 

multiplier asHociatcd with thrust direction, kg 

yaw attitude, rad 

multiplier associated with trajectory constraints 
dummy variable, seo 

adjoint multiplier associated with velocity, 
kg-seo 2 / nt 

adjoint multiplier associated with position, 
kg-scc/m 

adjoint multiplier >r ■■•tr.M* «.»d with mass, see 
adjoint multlplivt associated with t, kg 
declination or outgoing asymptote, rad 
pitch attitude, rati 
longitude of lino of nodes, rud 


3, Spurlock, O, F. andTcrcn, F,, "Optimum haunch 
Trajectories for Iho ATS- F Mission," AIAA Pn()or 
70-1031, Santa Barbara, Calif., 1070, 

<1. Teren, F. and Spurlock, O, I'., "Optimum Three- 
Dimensional Launch Vehicle Trajectories with Atti- 
tude and Attitude hate Constraints, " TN D-0U7, 
1000, NASA, 

Table I, - Baseline 'fug characteristics 


Tug initial weight (Including payload), kg . , . ,2b 022 


Tug burnout weight, kg 2707 

Engine thrust, N 723 

Engino specific Impulse, see , , , , ‘100 

Performance reserve, percent of 
specific impulse 2 


Table II, - Baseline mission characteristics 


Mission energy, (km/ecc) 2 ............. 12 

Unngo of declination, dog 0 to 31 

Paylond-soparnllon const time, min , , 10 

initial SS orbital allhudo, km . , 163 

Initial SS orbital inclination, deg ........ .28, C 


Subscripts: 


a npogco 

0 exhaust 

II hardware 

1 variable Index 

o processed SS orbit at time T^ 

p porlgco 

Superscripts: 

• time derivative 

vector 

' unit vector 
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Figure 1. - Pseudo right ascension of outgoing asymptote 
measured counterclockwise in equatorial plane from 
initial Space Shuttle orbit ascending node - without 
perigee propulsion. 
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Figure 2. - Argument of latitude measured In Initial 
Space Shuttle orbit - without perigee propulsion. 
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Figure 3. - Optimum Tug trip time from departure from 
Initial Space Shuttle orbit to return to Space Shuttle - 
without perigee propulsion. 
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Figure 4. - Nodal shift made during outbound burn and 
total nodal correction - without perigee propulsion. 
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Figure 5. - Inclination of hyperbolic or! 
sired outgoing asymptote - without pe 
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Figure 6. - Tug payload capability - without perigee 
propulsion. 
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Figure 7. - Velocity loss on outbound and inbound tra- 
jectory legs - without perigee propulsion. 



Figure 8. - Pseudo right ascension of outgoing asymptote 
measured counterclockwise in equatorial plane from 
initial Space Shuttle orbit ascending node - with perigee 
propulsion. 
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Figure 9. - Argument of latitude measured in Initial 
Space Shuttle orbit - with perigee propulsion. 
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Figure 10. - Optimum Tug trip time from departure from 
Initial Space Shuttle orbit to return to Space Shuttle - 
with perigee propulsion. 
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Figure 11. -Tug outbound first-burn time - with perigee 
propulsion. 
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Figure 12. - True anomaly at sftrt of second out- 
bound burn - with perigee propulsion. 
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Figure 13. - Nodal shift made during outbound burns and 
required total nodal correction - with perigee propulsion. 
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Figure 14, - Inclination of hyperbolic orbit containing de- 
sired outgoing asymptote - with perigee propulsion. 
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Figure 15. -Tug payload capability - with perigee propul- 
sion. 
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Figure 16. - Velocity loss on outbound and inbound tra- 
jectory legs - with perigee propulsion. 
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Figure 17. - Comparison of maximum Tug payload capa- 
bility, with and without perigee propulsion, with ideal 
performance. 
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